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The property of inhomogeneous turbulence in conducting fluids to expel large-scale magnetic fields in the direction of 
decreasing turbulence intensity is shown as important for the magnetic field dynamics near the base of a stellar convection 
zone. The downward diamagnetic pumping confines a fossil internal magnetic field in the radiative core so that the field 
geometry is appropriate for formation of the solar tachocline. For the stars of solar age, the diamagnetic confinement is 
efficient only if the ratio of turbulent magnetic diffusivity rjT of the convection zone to the (microscopic or turbulent) 
diffusivity r/in of the radiative interiour is rj-r/rii-n > 10''. Confinement in younger stars require larger r]T/riin- The 
observation of persistent magnetic structures on young solar-type stars can thus provide evidences for the nonexistence of 
tachoclines in stellar interiors and on the level of turbulence in radiative cores. 
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1 Introduction 

This paper concerns diamagnetism of stellar convective tur- 
bulence to show that the diamagnetic pumping can be im- 
portant for formation of the solar tachocline. 

Turbulent conducting fluids are known to expel large- 
scale magnetic fields in the direction of decreasing turbu- 
lent intensity. This effect of turbulent diamagnetism has 
been first found for inhomogeneous 2D turbulence (Zel- 
dovich [r957 ). For this case, the effective velocity of mag- 
netic field transport is proportional to the gradient of the tur- 
bulent magnetic diffusivity, Udin = — Vryx- The minus sign 
in the right shows the sense of turbulent magnetism: it is not 
para- but dia-magnetism so that magnetic fields are pushed 
away from the regions of relative high turbulent intensity. 
The diamagnetism is closely related to the magnetic field 
expulsion from regions of circular motion (Weiss il966l l. It 
was detected in 3D numerical simulations (Brandenburg et 
al. [19961 Tobiasetal.[l998l|200T; Dorch & Nordlund'200T] 
Thomas et al. 120021 Ziegler & Riidiger .2003) and labora- 
tory experiments with turbulent Uquid sodium (Spence et 
al. l2007i For 3D nearly isotropic turbulence the expression 
for the diamagnetic transport velocity becomes 



dia 



(1) 



(Krause & Radler |T980l l. The influence of rotation and mag- 
netic fields produce anisotropy and quenching of both turbu- 
lent pumping and diffusion (Kitchatinov |1988[ Kitchatinov 
& RUdiger 1992). In particular the strong magnetic quench- 
ing of the diamagnetism ~ B^^ for super-equipartition 
fields has been found. The isotropic parts of the field ad- 
vection and diffusion are still related by the Eq. ([l) when 
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Fig. 1 The radilal profile of the turbulent magnetic diffu- 
sivity in the solar convection zone after the model of Stix 
& Skaley (I1990l l. The diamagnetic pumping should be very 
strong near the base of the convection zone where the diffu- 
sivity almost jumps by orders of magnitude. 

rotation or strong magnetic field is imposed. Extensive nu- 
merical simulations - also as a test for the early theoretical 
SOCA-expressions - are due to Ossendrijver et al. ( 120021 ) 
and Kapyla et al. (|2006] |. 

The turbulent pumping can be significant in various as- 
trophysical contexts. Its strong effect is known in models of 
the galactic dynamo where the pumping may exceed the az- 
imuthal a— effect so that the dynamo can even switch off. In 
the present paper we suggest its importance also for the Sun 
and solar-type stars with external convection zones where it 
can participate in formation of the tachocline immediately 
below the base of the convection zone. 

Figure [T] shows the radial profile of the magnetic dif- 
fusivity in the solar convection zone estimated with the 
mixing-length relation rj^ — u'£/3. The diffusivity varies 
sharply near the base of convection zone. Its gradient repre- 
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sents a velocity of the downward diamagnetic pumping ([T]l. 
The amplitude of the velocity is up to 50 m/s exceeding both 
the amplitude of the a-effect (1-10 m/s) and the meridional 
flow in the convection zone ( < 10 m/s). 

The tachocline can be explained by an effect of a weak 
internal magnetic field of the solar radiative core (Riidiger 
& Kitchatinov 1 19971 120071 1 . The magnetic tachocline the- 
ory does not strongly restrict the magnetic field amplitude. 
Even a weak field well below 1 Gauss can produce the 
tachocline. The tachocline theory is, however, very sensitive 
to the field structure. The internal field should be almost to- 
tally confined inside the radiative core in order to produce 
the tachocline (MacGregor & Charbonneau l 19991 Kitchati- 
nov & Rudiger l2006l l. The downward diamagnetic pumping 
near the convection zone base can keep the internal mag- 
netic field confined in the radiative core thus providing the 
field geometry necessary for the tachocline formation (Ga- 
raud |2007l Garaud & Rogers |2007| |. The diamagnetic con- 
finement is estimated in this paper almost linear and in large 
scales. The main astrophysical question is how long the tur- 
bulent pumping would need to produce a confined magnetic 
field from an opened field configuration. Note that the dif- 
fusion time in the solar interior is longer by many order of 
magnitudes than the diffusion time in the convection zone. 
If the time to produce a confined magnetic geometry would, 
e.g., equal 10 Myr then - if the tachocline is magnetic by 
origin - the tachocline would only exist in older stars. On 
the other hand, the confinement mechanism via meridional 
flow with amplitudes of order 10 m/s works much faster 
(Kitchatinov & Riidiger [20061 ). There is thus the hope that 
future asteroseismology observations will decide which of 
both the mechanisms is really acting. 

Also the evolution of external fields in the convection 
zone under the presence of the diamagnetic effect must be 
reconsidered. This process may have consequences for the 
overshoot dynamo models (see Gilman 1992; Belvedere, 
Lanzafame & Proctor lT991t Riidiger & Brandenburg |1995l l 
although much stronger dynamo-fields may suppress turbu- 
lent advection of the field. 



2 Equations 

Consider the magnetic field in a sphere with an inner core 
where magnetic diffusivity rjin is relative low and an outer 
spherical shell with the diffusivity t^t ^ Vin- The diffu- 
sivity varies smoothly between these two values around the 
interface at r = 



1 , 

?7 = ?7in + 2 (^T - Vu 



1 + erf 



(2) 



where erf is the error function and defines the width of 
the transition layer The evolution of the large-scale mag- 
netic field is governed by the diffusion equation 



dB 

'dt 



which includes the diamagnetic transport with the effective 
velocity ([T]i- A vacuum boundary condition is applied at the 
top, r = R. 

Equation (O describes the decay of the field. We expect 
the decay being sufficiently slow for the internal fields to 
survive in the stellar radiation cores on time scales of Gyrs 
(Cowling |1945l l. With the diamagnetism included the mag- 
netic fields are expected to be confined in the core. A solu- 
tion of the eigenvalue problem 

B(r,i) = exp(cri)6(r), (4) 
provides both the decay time r — —a^^ and the field ge- 
ometry. The eigenmodes with the smallest decay rate is the 
most significant one. The degree of confinement of an ax- 
isymmetric poloidal field, 

' (5) 

r sm J 

can be described by the escape parameter 

H = rj, (6) 

(Kitchatinov & Riidiger 120061) . Here, r,9,(j) are standard 
spherical coordinates, and is the azimuthal unit vector. 
The parameter (|6]) measures the ratio of the characteristic 
values of magnetic flux through the interface to the flux in 
the core. The maximum value of (50 = 1 corresponds to 
an open field structure. The smaller 5(j) the more confined 
to the core the internal field is. A tachocline can be formed 
by the internal field if the escape parameter is sufficiently 
small, i.e. 

Scf> < 10-2, (7) 
(Rudiger & Kitchatinov l2007l l. 

The eigenvalue problem for the potential A of the 
poloidal field formulated in terms of the spherical harmon- 
ics, A ^ Ai {r)P^ (cos 9) sin 9, leads to the equation 
d / ^dAi\ + 



.A, 0^- 



A, 



(8) 



with the vacuum boundary condition 
at r = i?. 

We can solve Eq. ([8]) numerically for the continuous dif- 
fusivity profile (|2]i. But also an analytical solution can be 
found for a discontinuous profile of 77 (/id 0). In this case 
the solution of Eq. ([8]) for constant rj can be used, i.e. 

r \ , ^ , / r 



A, 



,(10) 



where is the Bessel function, r = —a^^ is the inverse 
eigenvalue, and a and b are free constants. Solutions with 
different sets of constants a and b apply to the inner core 
and the envelope. The boundary conditions 

" " " (11) 



[A] = 0, 



at r ~ Rn 



(3) 



can be applied at the interface. In the inner core, b must be 
zero to avoid a singularity at r = 0. The other constant a 
can freely be chosen, e.g. a = 1, to normalize the linear 
solution. Then, Eqs. (|9]l and (TT[ provide three conditions 
to define both the constants a and b in the outer shell and 
the eigenvalue. 
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3 Results and discussion 

3.1 Normal modes of the internal field 

Figure |2] shows the dependencies of the escape parameter 
dSJ and the decay time of the slowest decaying eigenmode 
of dipolar symmetry as a function of the diffusivity con- 
trast TjT I rjin between the inner core and outer envelope. The 
decay times are given separately for the analytical solution 
and numerical solutions for two values of transition width 
h. The trapping of the field in the radiative core grows with 
the diffusivity contrast. Figure [3] shows the structure of the 
eigenmodes. 
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Fig. 2 Dependence of the escape parameter Q and the 
decay time of the most long-living dipolar eigenmode on 
the diffusivity contrast ?7T/f?in- Full line corresponds to an- 
alytical solution for discontinuous change of rj. Dashed and 
dotted lines show the numerical results for smooth profiles 
© with /id = O.Oli? and = 0.02R respectively. All 
three lines overlap in the left panel. 

The tachocline can be formed by the internal field if the 
escape parameter is sufficiently small (|7]). After Fig.|2] this 
is the case if rjT/rjia > 10^. For an eddy diffusivity rj^ ^ 
10^'^ cm^s~^ in the convection zone, this inequality means 
that 

r?in < 10*^ cm^s^^ (12) 

in the tachocline. Hence, the tachocline can be only mildly 
turbulent if it is magnetic by origin. 

The internal field can also be confined by a global 
meridional flow penetrating the radiative core from con- 
vection zone (Kitchatinov & Riidiger |2006|. This type of 
confinement also can be efficient only if rjin is sufficiently 
low with the same upper bound (fT2l i as for the diamagnetic 
confinement. The reason for this coincidence might be that 
the confinements by meridional flow or inhomogeneous tur- 
bulence represent basically the same mechanism on differ- 
ent spatial scales. The condition ( fT2] l for the field expulsion 
from the region of motion to be efficient in the Sun does not, 
however, depend on the scale at all. 




Fig. 3 Poloidal field eigenmodes. The diffusivity contrast 
varies as riT/rjin = 10, 10'^, 10^ from the left to the right. 
The dotted circle shows the interface r ~ i?in. 



3.2 Internal field evolution towards a core-confined 
geometry 

A freely decaying field eventually approaches the most 
long-living eigenmode. To find the characteristic time scale 
for the field convergence to confined geometry of the eigen- 
modes the field was evolved in time with Eq. ^ starting 
from an initial field of open structure. Figure |4] shows the 
dependencies of the escape parameter (|6]l on time for sev- 
eral diffusivity contrasts. The S(j) decreases in time but even- 
tually approaches a constant escape parameter of the corre- 
sponding eigenmode. The confinement proceeds on the long 
time scale R^/rim of the internal diffusion despite the dia- 
magnetic pumping is fast. Why the confinement is so slow 
can be understood by observing the evolution of the mag- 
netic field structure. 
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Fig. 4 Escape parameter (|6]l as function of time for 
the runs starting from an open field structure. The lines 
are marked by the corresponding values of the diffusiv- 
ity contrast. The upper scale gives the physical time for 
the Sun computed for the microscopic diffusivity 771,1 = 
3 X 10^ cm^/s. 
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Figure |5] shows the initial phases of the confinement. 
After one (short) diffusion time, B? /rn, the strong con- 
centration of the field in the upper radiative core is already 
formed. The strong downward increase of the poloidal field 
across the interface is already in balance with diamagnetic 
pumping. The confinement slows down after this balance is 
achieved. The high concentration of the field immediately 
beneath the interface does not allow a high degree of con- 
finement even with a strong decrease of the field from the 
upper core to the lower convection zone. The field must be 
further smoothed over the entire core to attain the eigen- 
mode structure with high confinement. This smoothing is, 
of course, slow. Nevertheless, the internal fields can be con- 
fined already in young stars of the age of several tens of 
million years if the diffusivity contrast between the convec- 
tion zone and the radiative core is sufficiently high (Fig.|4|. 
Our calculations, therefore, predict that even rather young 
stars may akeady possess tachoclines. 

Observations can probe whether the solar-type stars of 
young clusters possess persistent magnetic structures which 
do not considerably change for years. If the fields rooted 
in the radiative core are not confined there, they will pene- 
trate to the surface with no remarkable variations. If persis- 
tent magnetic structures are found then the internal field is 
not confined and a tachocline cannot be formed. Simultane- 
ously, that would indicate a certain level of turbulence in the 
stellar cores. 
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